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This publication provides a coherent treatment for the reactor neutrino flux uncertainties suppres¬ 
sion, specially focussed on the latest 0 13 measurement. The treatment starts with single detector 
in single reactor site, most relevant for all reactor experiments beyond 9 13 . We demonstrate there 
is no trivial error cancellation, thus the flux systematic error can remain dominant even after the 
adoption of multi-detector configurations. However, three mechanisms for flux error suppression 
have been identified and calculated in the context of Double Chooz, Daya Bay and RENO sites. 

Our analysis computes the error suppression fraction using simplified scenarios to maximise relative 
comparison among experiments. We have validated the only mechanism exploited so far by exper¬ 
iments to improve the precision of the published 9 13 . The other two newly identified mechanisms 
could lead to total error flux cancellation under specific conditions and are expected to have major 
implications on the global #13 knowledge today. First, Double Chooz, in its final configuration, 
is the only experiment benefiting from a negligible reactor flux error due to a ~90% geometrical 
suppression. Second, Daya Bay and RENO could benefit from their partial geometrical cancellation, 
yielding a potential ~50% error suppression, thus significantly improving the global #13 precision 
today. And third, we illustrate the rationale behind further error suppression upon the exploitation 
of the inter-reactor error correlations, so far neglected. So, our publication is a key step forward in 
the context of high precision neutrino reactor experiments providing insight on the suppression of 
their intrinsic flux error uncertainty, thus affecting past and current experimental results, as well as 
the design of future experiments. 

PACS numbers: 


INTRODUCTION 

Reactor neutrinos have been used for fundamental re¬ 
search since the discovery of neutrinos. The last decade 
have witnessed a remarkable reduction of systematic 
error in reactor neutrino experiments, by about one 
of order of magnitude, imposed by the high precision 
needed to measure f?i 3 by Double Chooz [I] (DC), Daya 
Bay [2v3] (DB) and RENO @J, a milestone for the world 
strategy of neutrino flavour research. RENO has released 
several updates of the $13 analysis in conferences here dis¬ 
regarded until publications follow. The reactor measure¬ 
ments are consistent with all measurements sensitive to 
$13 obtained via other techniques providing a coher¬ 
ent di 3 perspective, as obtained by global fit analyses [Zl- 
19]. Since the reactor $13 experiments precision is unri¬ 
valled, they are expected to dominate the world knowl¬ 
edge on 0 i 3 , likely, for a few decades to go. Hence, reactor 
systematic dominates much of the 013 world knowledge, 
as experiments reach their final sensitivities. The mea¬ 
sured 0 i 3 (and its uncertainty) is expected to play a crit¬ 
ical role to constrain, or measure, still unknown neutrino 
oscillation observables, such as CP-violation and the at- 
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mospheric mass hierarchy M- 

To maximise the sensitivity to 0 i 3 , reactor experiments 
were forced to conceived experimental setups where flux, 
detection and background systematics are controlled to 
the unprecedented level of a few per-mille each contribu¬ 
tion (<1% total). The statical resolution is boosted by 
using multi-reactors sites. The unprecedented precision 
achieved is experimentally very challenging, therefore the 
redundancy 0 i 3 -experiments is critical, specially if their 
uncertainty budgets are complementary to offer maximal 
cross-validation. Despite some complementary, the reac¬ 
tor experimental setups are unavoidably similar and suf¬ 
fer from similar limitations, hence validation by different 
techniques would be important, although the precision 
needed is unattainable today. The fore-mentioned preci¬ 
sion improvement was obtained via multi-detector exper¬ 
imental setups, whereby, at least, two detectors are used 
for the reduction of the overall systematic budget since 
correlated systematics among detectors cancel out. This 
way, while the absolute systematics are the same, the 
relative systematics are much lower. The absolute sys¬ 
tematics are still dominant in any single-detector setup, 
such as DC (single detector), but also all past and, likely, 
most future reactor experiments. 

The systematics reduced by the multi-detector config¬ 
uration are: detection and flux systematics. Detection 
systematics benefit from dedicated detector design for 
them to be effectively identical , typically so, only upon 
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FIG. 1: The 2D geometry of the experimental setups of Double Chooz (left), RENO (middle) and Daya Bay (right) experiments 
is shown. The squares indicate the detectors and the circles indicate the reactor cores. The dotted lines depict the baselines 
between detectors and reactors, while distances are summarised in Tabic [TT] Notation on Daya Bay: ADI, AD3 and AD4 
correspond to the here called near-1, near-2 and far, respectively. 


full calibration, thus implying the same (or negligibly dif¬ 
ferent) responses and composition (cross-section, proton 
number, etc). Flux systematics benefit from the fact that 
the near detector(s) is located closer the reactor(s) such 
that the flux modulation originating from neutrino os¬ 
cillations, 0 i 3 in this case, is negligible (or very small) 
compared to the far detector(s) located further away. 
The far is placed at the expected maximal oscillation 
deficit driven by the Am 2 (atmospheric) constraint by 
MINOS [5] and T2K [B|. The suppression of the flux 
systematic is highly non-trivial, being the main subject 
of this publication. While having effectively identical 
detectors suffices to reduce detection systematics, from 
~2.0% [2] to ~0.2% [2], just having near and far detec¬ 
tors will not necessarily provide a full cancellation of flux 
systematics - unlike what was originally thought. There 
are three mechanisms leading to possible flux system¬ 
atic reduction to be elaborated in detail within the pa¬ 
per. First, multi-reactor uncorrelated uncertainties can 
benefit from having several identical reactors. Second, 
the near-far geometry of the experimental setup could 
enhance ability for the near become an effective perfect 
monitor to the far; i.e. cancelling fully the flux system¬ 
atic error. And, third, the nature of the reactor uncer¬ 
tainties; i.e. whether correlated or uncorrelated among 
reactors; might be exploited, as measured by multiple de¬ 
tectors. Reactor systematics suppress from single detec¬ 
tor scenario, typically ^3% [2] to <1% for multi-detector 
setups. DC achieves an impressive ~1.7% p] for a single¬ 
detector setup using the Bugey4 data for the mean cross- 
section per fission normalisation of the fission; i.e. as 
an effective normalisation-only near. The most accurate 
reactor flux anti-neutrino spectrum predictions mm 
rely on the ILL uranium and plutonium isotopes input 
data mini. However, the final flux systematics are a 
combination of the reactor and spectral systematic errors 
and depends on depends on each experiment configura¬ 


tion, since the evolution of the fission elements depends 
on the running configuration of each reactor. Therefore 
flux systematics are expected, with current knowledge, 
to be the dominant contribution for some of the reactor 
013 experiments, such as DB, itself leading the world 0i 3 
precision. 

This publication develops a framework for the non¬ 
trivial propagation of the reactor flux uncertainties and 
their suppression in the context of multi-detector and 
multi-reactor experimental setups, providing mechanism 
for the improvement of the global 0i3 precision. Most 
our discussion stays generic on flux systematic propa¬ 
gation; i.e. no need for the specific experiment error 
breakdown, allowing easy relative comparison across all 
experiments. Our study cases, however, inspired on the 
specific reactor- 0 i 3 experiments configurations for maxi¬ 
mal pertinence. The core of our calculations is analyti¬ 
cal, but a cross-check was implemented using a dedicated 
Monte-Carlo-based analysis, is also presented. The dis¬ 
cussion starts from the simplest single detector config¬ 
uration, then evolving towards a general formalism ap¬ 
plicable to any multi-detector and multi-reactor setups. 
The numbers linked to the specific experimental setups 
are however only guidelines as they are obtained in sim¬ 
plified scenarios, again, to maximise comparability. 


FLUX INDUCED UNCERTAINTIES 

So far, the reactor 0i 3 -experiments appear to treat the 
reactor systematics similarly. However, all collaborations 
are lacking dedicated publications on the topic, so ques¬ 
tions about the coherence across collaboration remain. 
Below, we will briefly summarise the reactor flux sys¬ 
tematics information, as provided by the different col¬ 
laborations. Most experiments rely on PWR reactors, 
hence the anti-neutrino flux is dominated by four iso- 
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TABLE I: Reactor-0i3 experiments breakdown of detector 
flux la systematic error Ob ;2 HI- The Daya Bay and RENO 
errors quoted here are correspond to their uncorrelated con¬ 
tributions, as they relied on multi-detector setups. Instead, 
the Double Chooz errors includes also the correlated contri¬ 
butions, as quoted from single-detector analysis. Therefore, 
the Double Chooz error is expected to be overestimated and 
will be revised by the collaboration in future publications. 
Thermal power (Pth), fission fractions (a/) and spent fuel are 
considered. 



Pth 

(%) 

O' f 

(%) 

Spent Fuel 
(%) 

Total 

(%) 

Double Chooz 

0.5 

0.9 

included 

1.0 

Daya Bay 

0.5 

0.6 

0.3 

0.8 

RENO 

0.5 

0.7 

unknown 

0.9 


topes: 235 U, 239 Pu, 241 Pu and 238 U (ordered by contribu¬ 
tion relevance). 

In general, flux systematics uncertainties are typically 
broken down into three terms: thermal power, fission 
fractions and spent fuel, as summarised in Table [I] The 
thermal power history is usually provided by the elec¬ 
tricity company exploiting the reactor. The associated 
error is related to the measurement method, the preci¬ 
sion of the installed sensors and the employed sensors 
calibration m • Each experiment has to precisely study 
the way how the thermal power is measured in order to 
estimate the corresponding error and the possible corre¬ 
lations among all the reactors involved, as typically, the 
same measurement techniques are applied to all identi¬ 
cal reactors in a power plant. The fractional fission rates 
have to be computed through reactor simulations, while 
their uncertainty estimation is not simple as they de¬ 
pend on the reactor model and the approximations used 
for each setup [16j. After each cycle, the reactors are 
stopped for refuelling, typically once per year in a few 
weeks lasting operation. In this operation, part of the 
fuel assemblies are exchanged with new ones. The spent 
fuel are stored in the dedicated pools, typically located 
next to the reactor sites with slightly different baselines. 
The spent fuel long lived fission products can generate 
a small fraction of anti-neutrinos above the inverse beta 
decay threshold m- Therefore, each experiment has to 
evaluate also the contribution from spent fuel and its un¬ 
certainty. 

DC, until now operating in single-detector mode (DC- 
I), quotes its reactor uncertainties as fully correlated be¬ 
tween reactors, which is the most conservative approach 
for this configuration as it will be demonstrated later on. 
The values for the thermal power and the fission fraction 
are 0.5% and 0.9%, respectively [lj. The spent fuel con¬ 
tribution was estimated within the dominant 1.7% flux 
systematic error quoted for all DC-I publications. The 
DC-I errors quoted account for both correlated and un¬ 
correlated contributions, thus they are overestimated rel¬ 


ative to its multi-detector configuration where only the 
uncorrelated errors are to be considered, as quoted by 
Daya Bay and RENO. DC is expected revise these num¬ 
bers in forthcoming the multi-detector DC (data taking 
started in 2014). However, our description handles er¬ 
ror suppression in relative, thus avoiding absolute error 
estimation and/or premature discussion, left for the pub¬ 
lications by the experiments. 

DB and RENO, instead, quote systematics for their 
multi-detector configurations; i.e. only the detector un¬ 
correlated components. DB quotes, as thermal power 
and fission fraction systematics, 0.5% and 0.6% respec¬ 
tively [2]. An extra contribution from the spent fuel is 
estimated to be 0.3%, increasing the overall systematic to 
0.8% (2]. DB presented a reduction of the total system¬ 
atic uncertainty to 0.04% on its ratio between observed 
over predicted rate upon optimisation, allowing modula¬ 
tion of the different near contributions relative to the far. 
This analysis, although not used for the measurement of 
013 , raises unsettled debate about the physical meaning 
of such a modulation, a priori fixed by the geometrical 
configuration of the site. For completeness, we have repli¬ 
cated such a result in this publication, however we disre¬ 
gard it fully for discussion on the precision of $ 13 . RENO, 
on the other hand, quotes 0.5% and 0.7% 0, respectively, 
for the thermal power and fission fractions, leading to an 
overall total reactor uncorrelated flux systematic uncer¬ 
tainty of 0.9%, in consistent agreement with DB. RENO 
quotes no error for the spent fuel, so it is neglected here¬ 
after, however this is somewhat unexpected since its con¬ 
tribution is expected to be non-negligible given the large 
number of reactors, like in the case of DB. 

THE ISO-FLUX CONDITION 

An important consideration for the reduction of reac¬ 
tor flux systematics is the geometry of the experimental 
setup, illustrated in Fig. [TJ which might lead to total er¬ 
ror cancellation. The location of the detector strongly 
depends on the overburden topology of the site, critical 
for background suppression. In a multi-detector exper¬ 
imental site configuration, the flux uncertainties would 
cancel entirely, regardless of the nature of the uncertain¬ 
ties, if the relative contribution made by each given reac¬ 
tor to the total detected antineutrino flux is the same for 
all the detectors. This is a simple acceptance condition, 
as it will be demonstrated later on. If met, the near be¬ 
comes an effective perfect monitor of the far, thus the flux 
systematic uncertainty is null. This condition is, in fact, 
trivially possible in case of isotropic sources, like reactors 
and, unfortunately, typically impractical in meson-decay 
neutrino induced beams (i.e. decay in flight), where labo¬ 
rious dedicated efforts, including dedicated experiments, 
are needed for an accurate spectral extrapolation from 
near to far, like in the case of MINOS and T2K. Com- 
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TABLE II: The distances in meters between detectors and reactors used for calculation, as illustrated in Fig [T] In the cases 
of DC and RENO, we inferred from [1, :4|, while for DB the information is clearly provided in [2]. 


Setup 

Ri 

R2 

R3 

R4 

As 

Re 

Double Chooz 

far 

near 

997.9 

351 

1114.6 

466 





Daya Bay 

far 

1920 

1894 

1533 

1534 

1551 

1525 

near-1 

362 

372 

903 

817 

1354 

1265 

near-2 

1332 

1358 

468 

490 

558 

499 

RENO 

far 

1556.5 

1456.2 

1395.9 

1381.3 

1413.8 

1490.1 

near 

667.9 

451.8 

304.8 

336.1 

513.9 

739.1 


plex multi-reactor sites experiments might not meet iso¬ 
flux condition. However, the iso-flux condition might be 
fulfilled partially, thus yielding partial suppression of sys- 
tematics. As the iso-flux condition is not fully fulfilled 
by any of the reactor-(?i 3 experiments, the estimation of 
the reactor systematics is not straightforward. So, it is 
mandatory to account for the differences among the reac¬ 
tor compositions fluxes and their corresponding system¬ 
atic errors. In fact, this will be investigated and quanti¬ 
tatively estimated in the following sections for each ex¬ 
periment, using some simple approximations. 

REACTOR UNCERTAINTY SUPPRESSION 

Let us consider a general experimental setup where one 
or more detectors D measure antineutrino fluxes gener¬ 
ated by N r reactors. For each detector, the incoming 
flux $£> represents a superposition of the fluxes emitted 
by Nr reactors having the magnitudes weighted by 
the corresponding solid angles subtended between the 
detector and each reactor. We express this as 

n r 

$£> = ^ ^DRi^Ri ( 1 ) 

Ri—Rl 

where index Ri runs over all the reactors (Rl, R2,... ,Nr) 
per site and the index D stands for the detector(s), where 
the near(s) and far(s) will be indicated later on by n and 
/, respectively. 

The accuracy and the precision of the geometrical solid 
angles is carefully controlled by the experiments by ded¬ 
icated surveys. The relevant distances are summarised 
in Table |TT| The associated uncertainty is generally so 
small relative to other terms, that we can neglect them 
for our analysis. We can also neglect the effect of neu¬ 
trino oscillation, as we consider the un-oscillated integral 
flux prediction. 

Before going through the specific cases, let us define 
three quantities to be used throughout the formalism. 


They characterise the systematic error associated to the 
emitted anti-neutrino flux contributing to the error of 
the measurement made by detectors. Sri represents the 
relative error of the emitted anti-neutrino flux per sin¬ 
gle reactor core Ri', i.e. the relative error of &Ri- Sd 
is the flux reactor uncertainty, as measured by the ex¬ 
periment; i.e. the relative error of <E>£>. For an experi¬ 
mentalist considering a precise measurement of $ 13 , the 
most important quantity is Sd, but the raw input is Sr. 
The ratio between Sd over Sr provides a measure for the 
effective error suppression fraction (SF) obtained upon 
error propagation, which is the core parameter used for 
error suppression characterisation. This factor reflects 
the ability of each experiment to reduce the overall re¬ 
actor uncertainty relative to the simplest case of one de¬ 
tector with one reactor configuration (i.e. Sr), where no 
cancellation is expected. The convention use is that the 
smaller the SF, the smaller the final reactor flux error 
systematic as measured at the detector (Sd), so SF is 
ranges within the interval [0,1]. The extreme values 0 
and 1 stand for total suppression (Sd = 0 , regardless of 
Sr) and no suppression (Sd = Sr), respectively. 


In this section, our formalism focuses on the propa¬ 
gation of the flux uncertainties, as characterised by the 
effective SF starting from Eq. |T] as applied to the spe¬ 
cific configuration of each setup and considering the same 
contribution for each individual reactor and ignoring the 
nature of each error component. Hence, most of our de¬ 
scription is done in relative; i.e. no necessity to introduce 
absolute error discussion; providing the ideal framework 
for highlighting the main features of each experiment 
in a comparable basis. As convention, we took the ac¬ 
tual experimental setups (next sub-sections) to illustrate 
our discussion. However, most of the discussion remains 
generic and figures will aim to be exemplify the general 
features. To simplify the description flow, some calcula¬ 
tions are only shown in the Appendix. 
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Single Detector: the DC-I Configuration 

This case is best illustrated by the DC-I (far only 
phase), where Nr stands for 2 reactors, Rl and R2. 
Propagating the errors on $£> (Eq. [I]) for the sole detec¬ 
tor involved, we obtain the corresponding SF expressed 
as 


SF 2 = 1 — 




( 2 ) 


{^Rl&Rl + £Ir2^R2) 2 

where k stands for the correlation factor between the gen¬ 
erated fluxes defined (see Eq. A34) as 


k = 


1 + 


s c B -si 



(3) 


where <5)) and 1 5r are, respectively, the correlated and the 
uncorrelated components of the total reactor flux uncer¬ 
tainty per individual core, defined as {Sr ) 2 = {Sr ) 2 + 
{Sr) 2 . Then, we can rewrite Eq. [ 2 ] as 

cp2 _ 1 + 

SF - (1+ !)*)" <4> 

using a term proportional to the difference between the 
reactor fluxes weighted by the corresponding solid angles 


Q _ ttR2$R2 ^ L 2 r1 ^ R2 
R “ L%^r 1 1 j 

where Lr\ and Lr 2 are the distances between the detec¬ 
tor and the reactors. The proportionality between the 
solid angles and the corresponding detector-reactor dis¬ 
tances is given by the inverse-square law. The results 
provided by Eq. [4] remain the same for a different choice 
of the reference reactor in Eq. [5j more precisely consid¬ 
ering the ratio {QriQri)/{£Lr 2 $R 2 ). Indeed, this case 
is equivalent to a transformation Q^r —> 1 /kldR which 
leaves Eq. [I] unchanged. 

The implications of Eq. [4] are best illustrated in Fig. [2] 
where the SF is shown against the reactor flux asym¬ 
metry between both Cliooz reactors, defined as ($^2 ~ 
and the reactor uncertainty type 
asymmetry, defined as ( S c — 8 U )/{S c + S u ), characterising 
the fraction of reactor error (un)correlation generically. 
There are two interesting limit cases in Fig. [2] where there 
is no error reduction (i.e. SF = 1): 

• when the errors are maximally correlated, repre¬ 
sented by the condition (<5 C — 8 U )/{5 C + S u ) —► 1.0 
(or k —> 1 in Eq. [ 4 ]) 



Flux Asymmetry 


FIG. 2: The suppression fraction (SF) for DC-I, coloured 
coded in the z-axis, is shown against the reactor fluxes asym¬ 
metry, defined as ($H 2 — < f > m)/( < f > .R 2 + $hi), and the uncer¬ 
tainty correlation asymmetry, defined as (<5 C — 5 U )/{5 C + 5 U ). 
No reduction of uncertainties is indicated by an SF=1.0 (red). 
Whenever one reactor dominates (i.e. the reactor flux asym¬ 
metry tends to ±1.0) or the reactor uncertainties are mostly 
correlated (i.e. (<5 C — S U )/(S C ± S u ) +1.0), there is small 
or even no SF. The maximal suppression (violet) is obtained 
whenever the reactor uncertainties are fully uncorrelated and 
both reactors are delivering similar fluxes. The expected SF 
minimum is l/y/2 (~0.7), only achievable for a slight offset 
to R2 to compensate the fact that the far and Rl are closer 
together, as shown in Fig. [l] 


• when either reactor is working, represented by the 
condition ($ fl2 - $ri)/($ri + $^ 2 ) -+ ±1.0 


this is because those cases cannot be effectively distin¬ 
guished from the case of single detector with a single re¬ 
actor scenario. If two reactors uncertainties are totally 
correlated, this is like if it was one effective reactor, re¬ 
gardless of the geometry of the setup. 

Conversely, some error reduction (SF <1) is expected 
elsewhere: whenever the source is consistent made up of 
two independent (i.e. uncorrelated errors; i.e. k = 0) 
reactors working at about the same contribution. Error 
suppression is expected to be maximal (i.e. SF minimal) 
when the two reactors have totally uncorrelated errors, 
since those reactor errors scale with l/\/2 (~0.7). The 
slight asymmetry in the reactor flux asymmetry (x-axis) 
is a consequence of the distance far-Rl being slightly 
closer than far-R2, as shown in Fig. [l] While rather 
trivial, the one detector case is important for many ex¬ 
periments and it is remarkably illustrative in the context 
of this discussion, since the fore-mentioned description 
inverts when considering a multi-detector sites, as de¬ 
scribed next. 
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Multi-Detector: the DC-II Configuration 

Let us consider now the SF behaviour for the DC-II 
scenario with two detectors with two reactors. In this 
case, it is more reasonable to characterised the SF asso¬ 
ciated to the near-far ratio , most relevant for oscillation 
analyses, as indicated by 


Ratio(//n) = (6) 

when C represents a constant such that the predicted 
un-oscillated flux at the far is obtained from the near 
(i.e. $ r j. red = C x <f> n ). This constant accounts for the 
difference in the baselines between the near and far sites, 
due to the isotropic emission of neutrinos in reactor. The 
near and far are respectively indicated by the subindexes 
n and /. Let us reconsider Eq. [l] now generalised to 
consider the contribution of the two detectors to compute 
the ratio $//$„. As in the case of the single detector 
case, we are interested in the expression for SF but now 
relative to the ratio between detectors, indicated by S n f, 
to the relative error of the flux emitted by each reactor, 
thus 


SF 2 = 



fv) 3 +fv') S - 2E S# 1 

\5rJ \5rJ SjiMf 


(7) 


where S n and S / are the relative errors of the fluxes reach¬ 
ing the individual detectors, respectively for the near and 
far. The first two terms in Eq. [7] are the same as the ones 
obtained for single detector configuration (Eq. [4]) . The 
new term corresponds to the covariant term computed as 


COv(<I> n , $/) = SxlttnRittfRi&ft! 

+ QnR2^fR2®R2 
+ k$Ri$R2(fl nRl^fR2 
+ ^nR2^fRl)) (8) 

where k has the same definition as in Eq. [3] for the error 
correlation factor. We introduce fl n R and ftfR parame¬ 
ters, as in Eq. [5j to obtain now 

cov($n,$/) = My ^nR^fR + k ($lnR + &f r) /qn 

(1 + f2/jj)(l + flnfl) 

such that the previous Eq. [7] now becomes 

qp2 _ 1 + ^nR + 2fcf l n R 

(i + n nR ) 2 

1 T r T 2k f r, 

+ (i + n fR ) 2 

0 1 + ^nR^fR + k(fl n R + f IfR,) . 

(i + n f R)(i + n nR ) 1 j 


representing SF general expression for the near-far setup 
with two reactors. 

As before, fl n R and QfR are independent from any 
factor which multiplies both solid angles or both reac¬ 
tor rates, thus the weight parameter from Eq. [6] does 
not contribute to the expression from Eq. [To] This 
expression is also independent from the transformation 
floR —> 1 /^dr (with D = n,f) given by a different way 
of defining the fraction terms in Eq. [5j By definition, the 
iso-flux condition implies that the relative antineutrino 
fluxes from reactors is the same for both detectors. This 
can be mathematically represented as 

ttnR = ftfR (11) 

The implication of Eq. [lOjare best illustrated in Fig. [3] 
where the evolution of the SF is shown. Fig. [3]-Left plot 
shows SF evolution relative to both the reactor accep¬ 
tance asymmetry, defined as (flR 2 — Qri)/(Qr 2 + Qri), 
and the reactor uncertainty type asymmetry, defined as 
(5 C — 6 U )/(5 C + S u ). The iso-flux condition is met when¬ 
ever the reactor acceptance asymmetry is ~0, leading to 
a SF to fully cancel, regardless of the nature of reactor 
uncertainties, demonstrating so the near is a geometrical 
perfect monitor to the far. This is the same as impos¬ 
ing the condition in Eq. [Til in Eq. [lOj The DC geometry 
is such that the iso-flux-ness is partially met along the 
projection (I 1r 2 — SIr\)/{£Ir 2 + Hfli) ~ —0.28, hence a 
major reduction of SF is obtained where SF can, at most, 
be ~0.12; i.e. M0% of the original uncertainty is thus 
suppressed. Fig. §R ight plot shows SF evolution rela¬ 
tive to both the reactor power flux asymmetry, defined 
as (d>i? 2 — < hfli)/(^fl 2 + < i ) i?.i), and the reactor uncertainty 
type asymmetry, as defined before. The effect of the reac¬ 
tor flux asymmetry can only exemplified in the case of DC 
with two reactors. Total cancellation of SF also is found 
if the uncertainties are reactor correlated maximally (i.e. 
k = 1) or either reactor is off. This is expected because 
both such conditions are equivalent as having one effec¬ 
tive reactor as source, thus perfectly monitored by the 
near, regardless of the experimental setup geometry or 
the type of the uncertainty. As expected, this conclusion 
is general and true independently of how many reactors 
are considered. This implies that, for example, whenever 
DC runs with 1 reactor and both detectors, the flux er¬ 
ror is zero for that data set. It is worth noting that the 
pattern shown for DC-II, illustrated in Fig. [3]- Right, is 
exactly the opposite to the one exhibited by DC-I, illus¬ 
trated in Fig. [2j where no suppression (SF=1) is obtained 
in the case of total correlated reactor uncertainties. In 
brief, the DC-II performance, SF ~0.12, exemplifies the 
best multi-detector error cancellation due its almost iso¬ 
flux geometry where almost ~90% of the SF is due ge¬ 
ometry. The remaining SF fraction (~30%) is caused by 
the 1/y/N r scaling of the remaining totally uncorrelated 
error, where Nr = 2. 
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Acceptance Asymmetry 


Flux Asymmetry 


FIG. 3: The suppression fraction (SF) for DC-II is illustrated coloured coded in the z-axis. Left shows the evolution of SF 
against both the reactor acceptance asymmetry (x-axis), defined as ( 0 _H 2 — Qri)/(Qr 2 + flfii), and the reactor uncertainty 
type asymmetry (y-axis), defined as ( 5 C — 5 U )/(8 C + 8 U ). The iso-flux condition is met whenever the acceptance asymmetry 
is ~0, leading to a SF —>0 (violet), regardless of the nature of reactor uncertainties, thus demonstrating that the near is a 
geometrical perfect monitor of the far. The DC geometry is such that the iso-flux-ness is partially met along the projection 
( Qr 2 — Diji)/(f2ij2 + Dfli) ~ —0.28, hence SF is, at most, ~0.12. Right shows the evolution of SF against both the reactor 
power flux asymmetry (x-axis), defined as (<&R 2 — &ri)/(&r 2 + "Fm) (he. the flux difference between reactor R1 and R2), and 
the reactor uncertainty type asymmetry (y-axis), as defined before. Note that the maximal value possible of SF is the same 
~0.12, as explained. Total cancellation (i.e. lowest SF) is found if the uncertainties are reactor correlated maximally or either 
reactor is off. This is because both cases are mathematically identical as to having one effective reactor as source, perfectly 
monitored by the near, regardless of the experimental setup geometry. 


Multi-Detector: the RENO Configuration 


Generalising Eq. [To] to account for a larger number of 
reactors (Nr = 6), RENO is becomes our next case with 
two detector sites and Nr = 6 , as illustrated in Fig. |T| 
The new general expression is 

CI? 2 ( 1 + Sd 2 + 2 k(Sd + S' d ) ^ 

SF = M-- ) 

d=n,f 

o l + S nf + k(S n + S f + S ' nf ) /ioN 
(1 + S n )(l + Sf) 1 j 

where we have introduced the following notation 


Nr -i 


q' _ oW qC?) . 

and £l d R is a generalization from the Eq. [HJ to obtain 


,(* 

dR 


^dRi^Ri 

^dRl^Rl 


LdRl^Ri . 


1 = 2...Nr-1 (13) 


where reactor R1 has been arbitrary chosen as reference 
reactor. Nevertheless, choosing a different reactor leave 
Eq. [12] unchanged. 

The iso-flux condition, originally in Eq. El is now more 
generally expressed as 


n r -i 


Sd= E 


i> 1 
Nr- 1 


^ = E 

i> 1 

Nr- 1 

s'd= E 

i,j>l m ,i^j 

Nr- 1 

c „ _ \ o(») rd®) • 

~ 1 ^f R 1 

*>1 


tonR = to?R (14) 

which, like in the case of DC-II, would cause SF to go 
to zero, regardless of any other error dependence, upon 
imposition, if the geometry of the experimental site con¬ 
sidered allows. 

Upon introducing the geometry of RENO, we can sum¬ 
marise our finding as follows, 

• SF will obtained full cancellation if reactor errors 
are fully correlated, as demonstrated and illus¬ 
trated for DC-II (Fig. [3]), regardless of the experi¬ 
mental site geometry. 












• SF will be largest in case of the fully uncorrelated 
reactor systematics errors when they deliver sim¬ 
ilar fluxes, also as obtained for DC-II. However, 
an overall SF follows 1 /y/N R reduction, hence, SF 
will benefit sites with many independent (i.e. un¬ 
correlated errors) identical reactors. In a sites with 
different reactor types, the effective SF might be 
deteriorated. 

Mathematically, this can be seen when considering 
k = 0 and Nr —> oo into Eq. [12] The factors on 
the denominator, (1 + Sd) 2 and (1 + S n )(l + Sf) 
will increase faster than the corresponding factors 
at numerator, S,p respectively S n f , hence leading 
to cancellation. 

• RENO has some SF reduction due to some degree 
of iso-flux-ness. The amount of SF reduction is lim¬ 
ited, as expected, since the geometry of the site is 
not optimal to ensure that the near and far have 
a similar contribution across all reactors equally. 
The near sees, respectively, ~78%, ~16% and ~ 6 % 
from the different reactors pairs (from the nearest 
to the farthest), while the far sees a more even con¬ 
tribution from all reactors. Thus, the RENO site 
is not expected to meet the iso-flux condition fully. 
The total SF measured was estimated to be ~0.23 
in the full power scenario. 

• There is a small difference between the nominal 
power of two RENO reactors [4]. In our analysis, 
we took this into account by weighting the corre¬ 
spondent terms. 


explicit contribution of two near sites in the denominator. 
Hence, the new expression becomes 


Ratio(//n) 


/3$nl + 7 $n2 


(15) 


where /3 and 7 are two constants, whose values are 
fixed by the experimental geometry and the fluxes of the 
running reactors. If optimisation of the site, upon design, 
wanted to be considered, the effective SF of DB would 
have depended on them. For the prediction, the new flux 
is thus constructed from the combination of the two near 
sites measurements; i.e. $ n 2 - However, the near 

sites do not only monitor their respective reactors, but 
instead, they also see a fraction of the other reactors (not 
supposed to be monitoring); i.e. near -2 also sees ~6.5% 
of reactors R1 and R2 and near-1 see about, ~17% of the 
remaining reactors. This makes this ratio very delicate, 
since some degree of double-counting is unavoidable by 
the nears is unavoidable, while this is not present in the 
far. This will translate into a loss of the iso-flux-ness con¬ 
dition. The same conclusions listed for RENO are valid 
for DB (previous section). The overall SF estimated for 
this site is 0.18 benefiting for slightly better partial iso- 
flux-ness relative to RENO, despite the visual appealing 
geometrical symmetry of the RENO site. In the case of 
both RENO and DB, most of the SF (~60%) arises from 
the 1 /V Nr scaling due to their large number of identical 
reactors. 


Suppression Fraction Estimation via Monte-Carlo 


Unlike so far assumed, commercial reactors do not de¬ 
liver their flux constantly over time. Some unavoidable 
variations expected are due to refuelling, a total reactor 
stop once a year over a few weeks, and reactor burn-up 
effects, typically a few percent decrease in flux over the 
entire year. In addition, reactors might be run differ¬ 
ently every year, so each reactor cycle history is a priori 
expected to be unique. When that happens in site with 
many reactors, the effective SF can only deteriorate as 
the two above conditions will be varied away from its op¬ 
timal spot. This effect will be further studied and quan¬ 
tified later on. 

Multi-Detector: the Daya Bay Configuration 

Let us now consider the setup of the DB experiment 
now. As illustrated in Fig. [T] this site differs from the 
other on having 2 near sites to monitor the two sets of 
reactors geometrically grouped of the power plant. The 
near-1 monitors mainly reactors R1 and R2, while and 
the near -2 monitors all others. 

Using similar formalism till now, the only modification 
here is to build the Ratio(//n), used in Eq. | 6 j with the 


Together with the analytical formalism presented, we 
have analysed the reactor flux systematics using a Monte- 
Carlo based flux simulation method. Likewise, the full 
geometry and ingredients to each experimental sites are 
simulated. Each reactor is described by its thermal power 
(Pth) and fractional fission rates (a*, with k = 235 U, 238 U, 
23ri Pu and 241 Pu, the main core isotopes), while each de¬ 
tector is described by the number of free protons in the 
target and by the detection efficiency. For simplicity, the 
same mean au ’s have been used for all experiments, how¬ 
ever this is not totally correct, as different experiments 
not only have different reactors but also the fuel com¬ 
position, namely the enrichment of 235 U, is expected to 
be somewhat different. Since the detector correlated er¬ 
rors are canceled in multi-detector setups, we have prop¬ 
agated only the detector uncorrelated errors. 

The antineutrino flux emitted by a given reactor is a 
function of thermal power and the fission rates define 

$ = <f>(P th ,a k ). (16) 

As the error propagation formalism does not depend on 
the specific form of the function <F, we do not describe 
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it in detail. For a given binning of the antineutrino en¬ 
ergy spectrum, the error of the flux is expressed as the 
covariance matrix 

(17) 

where i and j stand for the energy bins, and M Pth = 

{Sp <5p ) and M°j = (6 l a 6 J a ) are the uncertainty con¬ 
tributions from the thermal power and the fission rates, 
respectively. As the error on does not depend on 
the antineutrino energy, M p,h is a diagonal matrix with 
6 l P =6 J P =S Pth . 

For a given reactor r, the error on <f> due to the thermal 
power uncertainty Sp th is computed as 

3p th = &r{P[h + a P th i a k) — ®r(Pth><*k), ( 18 ) 

where Op is the 1 -a error of the thermal power. 

A7“ is computed by propagating the covariance matrix 
of the fission rates, C ^ = pu^ak^on, being p the correla¬ 
tion matrix and er Q the 1-a error of the fission rate 


M% r = Y, D * C ki D Ji> (! 9 ) 

k,i 

where D stands for the derivatives defined as 


SF= (23) 
s f 

where 8* NF is the error on the flux at the far detector 
when computed as ratio(/n) x $ n . Such a SF is inde¬ 
pendent of the parameterisation of <f> and the specific 
values of crp th and cr a . 

In this calculation, the SF accounts only for the error 
suppression provided by the iso-flux-ness of the experi¬ 
mental setup. In order to be compared with the analyti¬ 
cal estimation, developed in the previous sections based 
on the S*/5p ratio, one needs to incorporate the uncor¬ 
related error reduction due to the number of reactors. If 
detectors fulfilled the iso-flux condition, such a reduction 
would be \/Np. As this is not the case for any of the 
experiments considered, an effective number of reactors 
(Npf) has to be computed. According to our simula¬ 
tion, the Np f values are 1.98, 5.95 and 5.78, for DC, 
DB and RENO, respectively. 

The SF obtained via Monte-Carlo, once corrected by 
the -Mp , are in remarkably good agreement with the 
ones from the analytic calculations: identical numbers 
reproduced to the 3 rd digit. This provides a mutual val¬ 
idation of the procedures presented. 

Reactor Time Variations: Refuelling Scenario 


D = d»;(q fc + tJfc, P th ) - (a fc , P th ) 


( 20 ) 


The a-related error on the integrated flux from a reac¬ 
tor R can be expressed as 


^ = ( 21 ) 

where £) stands for the energy phase space of bin i. 

For an experimental setup with Nr reactors, the final 
integrated error on the antineutrino flux is derived as 
the quadratic sum of the contributions from individual 
reactors, provided by Eq. [18] and Eq. [21] 


A* = 


\ 


Nr 




( 22 ) 


Using the procedure described above, we computed the 
central values ($" and $•') and uncertainties (S n and 5?) 
of the fluxes at both near (n) and far (/) detectors of DC- 
II, DB and RENO experiments. We also computed the 
far-to-near flux ratio, ratio (fn) = along with its 

error and, finaly, derived the error SF as 


Considering time variation of the flux delivered per re¬ 
actor will modify the effective SF behaviour described so 
far, as the effective site reactor power configuration de¬ 
parts from the full-power over-simplified scenario consid¬ 
ered so far. In this section, we shall consider the impact 
to the effective SF for each experiments upon reactor 
refuelling estimated analytically, which a more realistic 
experimental scenario since this is, by far, the largest 
expected reactor-to-reactor variation that can be mod¬ 
elled similarly to all sites considered. The effects asso¬ 
ciated to burn-up, fuel composition, variations caused 
by running operation constraints (unicity reactor cycle) 
will all be neglected, as they are expected to be both 
smaller in magnitude and site dependent, thus out of the 
scope of the generic description considered here. How¬ 
ever, more accurate estimations expected by each specific 
experiment should consider its impact carefully. There¬ 
fore, for simplicity, let us consider the following refuelling 
scenario: for a total period of one year, each reactor is 
stopped for two months for refuelling, while during this 
period the other reactors are running. In this scenario, 
we consider also fully uncorrelated reactors errors, rep¬ 
resenting the worst case scenario. The impact on SF 
due to one reactor going down is, obviously, site depen¬ 
dent. For DC-II, the overall SF is zero whenever one 
reactor is monitored by two detectors, regardless of the 
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TABLE III: SF values for different RENO reactor configura¬ 
tions: all ON and each off. Notation as shown in Fig. [l] 


ON 

R1 off 

R2 OFF 

R3 off 

R4 OFF 

R5 off 

R6 off 

0.227 

0.226 

0.265 

0.234 

0.285 

0.244 

0.206 


geometry, as demonstrated before. For RENO and DB, 
contrary to DC-II, we expect an effective deterioration of 
SF since both the effective iso-flux-ness condition (power 
symmetry geometry) and the number of reactors running 
worsens upon refuelling. 

Our results can be summarised as follows: 

DC-II: the year average SF is 0.08 (refuelling scenario), 
to be compared with 0.12 (full power). Any one 
reactor off periods will only benefit DC-II as it has 
null flux uncertainty, therefore the longer the better 
for systematics, although that will imply a loss in 
luminosity. 

DB: the tear average SF is 0.20, against the 0.18 ob¬ 
tained at full power. 

RENO: the results are shown in Table ITTT1 for each reac¬ 
tor, however the year average SF was found to be 
0.24, to be compared when considering all reactors 
(0.23). 

Finally, the refuelling scenario here described (one re¬ 
actor stopped at the time) is specially naive in the case of 
RENO and DB, as their data shows simultaneous stop of 
up 3 reactors for some time, hence the site power might 
swing up to ~50%. This will have implications, typi¬ 
cally, towards the deterioration of the effective SF to be 
computed for those sites, when integrating over all those 
effects, including fuel burn-up. That level of accuracy is 
beyond the scope of this estimation. 


Daya Bay Near Detectors Optimisation 

As indicated in Eq. |15[ combining flux measurements 
from different nears with different weights seems an ap¬ 
pealing consequence for the case of DB. This represents a 
fundamental difference with respect to the experimental 
configuration with a single near. Any scaling of a flux 
measurement in Eq. |T] is equivalent to a change to the 
baseline of the detectors. So, a difference between the 
values of /3 and 7 , in Eq. [l5j is equivalent to a change be¬ 
tween the relative baselines of the near sites. The authors 
think that such transformation might become problem¬ 
atic given two possible effects: 

• the near sites do not only monitor the closest reac¬ 
tors, but there is a sizeable fraction of flux coming 
from the further reactors. This implies that the 



FIG. 4: The SF variation of DB for the near sites relative 
weighting is illustrated. The minimum obtained is in agree¬ 
ment with DB published results. The case for all reactors at 
full power is indicated by the dashed line, while the maxi¬ 
mal variations obtained due to refuelling scenario (simplified 
case described in text) is indicated by the shaded region. The 
physical interpretation of this analysis and its application for 
meaningful reactor flux error suppression is a debate with no 
consensus within the field. Our result is here illustrated as 
validation of published results, but it will disregarded for fur¬ 
ther discussion. 


relative normalisation of each near are not fully in¬ 
dependent. 

• the effect of the neutrino oscillations measured by 
the near detectors are not fully negligible. 

If we ignored both physical conditions, we could con¬ 
sider the optimisation of the contribution of the near 
sites, as done by DB in Jj. As we have shown previ¬ 
ously, our formalism is transparent to any factor which 
multiplies all solid angles or all fluxes provided by re¬ 
actors with the same factor for a given detector. Thus, 
we consider only the relative difference of the weighting 
parameters as 


Nr 

4> n = <F P f rediction = $ ni + p$„ 2 = n nR i(p)<$>Ri (24) 

i 

where p = 7 //3 is the relative weight parameter and 


^nRi(p) — n\Ri T P^n^Ri (25) 

is a linear combination of the omega terms as described 


by Eq. 13 


Since both 4>„ and $/ are given by Eq. [TJ we can use 
directly Eq. [12] in order to get the expression for the 
corresponding SF. Inserting the values of the geometri¬ 
cal parameters, we show the variation of the SF relative 
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TABLE IV: The analytical estimated suppression fraction (SF) per experiment are here summarised. The SF estimated via MC 
are not shown, as they are numerical identical up to the third digit. SF(total) is shown in both the full-power (over-simplified) 
and the refuelling scenarios. The SF(total) is defined as SF(total) = SF(iso-flux) x SF (Nr) x SF (correlation), whose value 
is within the interval [0,1], implying total and no suppression, respectively. The potential exploitation of SF(correlation) 
suppression needs dedicated analysis beyond the scope of this publication. Hence the remaining error is assumed to be totally 
reactor uncorrelated, implying a SF (correlation)=1 (i.e. no suppression), being the most conservative scenario for multi¬ 
detector and multi-reactor experiments. Note, however, that this implies maximal suppression due to SF (Nr). The impact of 
SF(correlation), while not estimated for any specific case (being strongly reactor dependent) is generically illustrated in Fig. [5] 


Experiment 

Nr (via MC) 

SF (full-power) 

SF (refuelling) 

SF (iso-flux) 

SF (Nr) 

Day a Bay 

6 (~6.0) 

0.18 

0.20 

0.49 

0.41 

Double Chooz 

2 (~2.0) 

0.12 

0.08 

0.11 

0.71 

RENO 

6 (~5.8) 

0.23 

0.24 

0.59 

0.41 


weighting parameter in Fig. [4] for fully uncorrelated er¬ 
rors of the reactor fluxes. The dotted curve represents 
the variation of the SF when all reactors are operational 
and the continuous curve represents the mean value be¬ 
tween the values obtained when one reactor is off upon 
reactor refuelling scenario considered in the last section. 
The variations in such refuelling scenario are indicated 
by the grey area in Fig. [4] accounting for the different re¬ 
actors. There are two interesting cases to be highlighted 

p — 1: this point represents no relative weighting ap¬ 
plied on the near sites, for which SF is ~0.20. The 
spread of SF when considering the reactor refuelling 
scenario is ~ 2 %. 

p at minimum: the effective SF obtained is about 
~0.05. The spread due to refuelling is about ~ 6 %. 

Our result is consistent with the official result of 
DB J|, where they obtained /3 = 0.04 and 7 = 0.3 for 
the analysed data period, considering that one of the near 
contains twice the events than the other, we obtained p 
to be 3.38. Using MC, we also obtained respectively, 
0.04 and 0.3. Our calculation here, thus serves as a repli¬ 
cation. However, due to the difficulties on the physical 
interpretation of this, the corresponding SF is ignored 
any further in this publication. This is consistent with 
the fact that DB discards it for the measurement of 6 * 13 . 


SUMMARY & DISCUSSION 

We have identified, studied and quantified three mech¬ 
anisms inducing reactor flux uncertainty suppression 
in the context of multi-detector experiments in multi¬ 
reactor sites. We have quantified the integral error 
suppression by the suppression fraction (SF) analyt¬ 
ically (cross-checked via MC) using simplified exper¬ 
imental scenarios to allow coherent relative compari¬ 
son across experiments. SF can take values within 
[ 0 , 1 ], where the extreme cases stand, respectively, for 
total suppression (SF=0) and no suppression(SF=l). 
The three mechanisms can be characterised by their 


respective SF terms, since SF(total) is defined as 
SF (total) = SF (Nr) x SF(iso-ffux) x SF (correlation), 
where i) SF (Nr) is linked to 1 /VNr scaling of 
the remaining uncorrelated reactor error, ii) SF(iso- 
flux) is linked to the site iso-flux condition and in) 
SF (correlation) is linked to the nature of reactor errors. 
Both SF(iso-flux) and SF(correlation) could lead to total 
error suppression under specific site conditions. If those 
terms are to be exploited, those conditions must be care¬ 
fully evaluated and demonstrated by each experiments 
accounting accurately for all pertinent effects, although 
no experiment have ever done this. Our final results rely¬ 
ing simplified refuelling scenarios are summarised in Ta¬ 
ble IV This results are not expected to be used by exper¬ 
iments, but as mere guidelines for more accurate estima¬ 
tions to follow up. The contribution of SF(correlation) 
was not quantified for any specific experimental setup, as 
it deserves more careful treatment discussed below. 

The SF (Nr) has been actively exploited in publica¬ 
tions by experiments to improve our knowledge on $13 
under the assumption that the remaining error is is fully 
reactor uncorrelated. Typically, Nr refers to the number 
of effectively identical reactors per site. In the case of 
DB and RENO, this term amounts to ^60% suppression 
(6 reactors), whereas this yields only ~30% suppression 
for DC-II. Typically, SF (Nr) is propagated into the $13 
precision as a byproduct of the y 2 minimisation. 

The estimation of SF(iso-flux) has not been applied 
by any experiment so far. It is likely unpractically to 
be implemented via the y 2 minimisation formulation, in¬ 
stead calculations might follows the prescription here pre¬ 
sented. Once estimated, as demonstrated in this publi¬ 
cation, the SF (iso-flux) is expected to improve the so far 
published results by DB and RENO, providing an extra 
flux error reduction by up to ~50% and ~40%, respec¬ 
tively on 013. Due to the simplified conditions assumed 
for our calculations, our SF’s are expected to be slightly 
optimistic relative to those to be obtained by dedicated 
analyses by DB and RENO eventually. In the case DC- 
II, the SF(iso-flux) term is expected to yield a dramatic 
~90% error reduction since the iso-flux condition is al¬ 
most met. This makes DC-II the only $13 experiment 
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likely to benefit from a negligible flux error, compared 
to other systematics. DC has officially prospected a con¬ 
servative 0.1% as flux error for DC-II [3], well within 
the analysis here presented. Therefore, DC-II final 6*13 
sensitivity is expected to be dominated its challenging 
background systematic,thus in maximal complementary 
to DB error budget. Since the flux error is dominant 
for DB, the hereby presented error reduction by ^50% 
translates into a significant improvement of the world 0 13 
precision by means of DB alone, but also via the envis¬ 
aged combination by all reactor experiments. 



FIG. 5: The variation of SF to the reactor uncertainty type 
asymmetry, defined as (5 C — 5 U )/(S C + 5 U ), is shown for DC-II, 
RENO and DB. All experiments have assumed, so far, their 
errors to be reactor uncorrelated; i.e. (<5 C — 5 U )/{5 C + 8 U ) = 
—1.0. DC-II benefits the best SF (~0.1) due to its almost 
iso-flux site. RENO and DB benefit mainly from the large 
number of reactor error suppression, but they also have some 
partial iso-flux matching hence benefiting from an extra up 
to ~50 error suppression so far neglected in 9 13 publications. 
Both full-reactor power (dashed lines) and a simplified re¬ 
fuelling scenario (solid lines) are shown, while the latter is 
expected to be a more accurate description of the reality. In 
the refuelling scenario, DC-II benefits from the null reactor 
flux error whenever only one reactor is running, while RENO 
and DB have the expected opposite trend. 

The SF(correlation) term has the potential to pro¬ 
vide further flux systematic error suppression, regard¬ 
less of the reactor site geometry, in the (unexpected) 
limit of full correlation across reactors, SF—>0 since the 
SF (correlation) term cancels, as illustrated in Fig. [ 5 ] As¬ 
sessing the correlation among reactor errors is a diffi¬ 
cult subject, having today no specific common prescrip¬ 
tion or consensus. Thus the only existing consensus is 
to adopt the most conservative scenario, implying two 
distinct cases: 

Single-Detector Case: maximal SF(total) is obtained 
if total correlation of reactor errors is assumed, 
as shown in Fig. [2] This is because both 


SF(correlation) and SF(N R ) terms provide no error 
suppression; i.e. SF=1 each. This the scenario is 
assumed for all DC-I publications and all single de¬ 
tector experiments, unless otherwise proved, thus 
affecting all past and and most future experiments. 

Multi-Detector Case: maximal SF(total) is obtained 
if total uncorrelation of reactor errors is as¬ 
sumed, as shown in Fig. §R ight. This implies 
SF (correlation)=1 and SF(N R )=l/y/N R , hence 
benefiting from maximal SF(N R ) reduction. This 
is the scenario expected to be assumed by DB, DC- 
II and RENO, until otherwise proved. 

Of course, reactor errors are unlikely to be neither to¬ 
tally correlated or totally uncorrelated. The fact that 
only those extreme cases are considered in the literate is a 
mere demonstration of the lack of knowledge for a better 
handling. Beyond the current debate among experts on 
the subject, the promising exploitation of SF(correlation) 
will require strong reactor-type dependences to be ac¬ 
counted and justified carefully by each experiment in 
dedicated publications. This means that each experiment 
will have to analyse their respective reactors, thus pro¬ 
viding insight evidence of the error correlation behaviour. 
The thermal power contribution, typically indicated by 
Pth, depends mainly on the uncertainty analysis of the 
internal reactor instrumentation data used for power es¬ 
timation, as provided by the reactor running company. 
Instead, the fission fraction evolution, indicated by a/, 
depends mainly on the simulation uncertainties analy¬ 
sis, including the assumptions and input parameters (fuel 
configuration, etc) used for the modelling and time evolu¬ 
tion. Both the instrumentation and simulations are very 
specific to each reactor type and, therefore, to the each 
experiment. Therefore, dedicated analyses are needed 
by each experiment to justify the delicate exploitation 
of SF (correlation), in the same way that this publication 
aims to illustrates the realisation of the unprecedented 
SF (iso-flux) exploitation. 

As a consequence of the unsettled complications be¬ 
hind the assessment of SF (correlation) for each experi¬ 
ment, our description here remains generic in its descrip¬ 
tion, as shown in Fig. [5] nonetheless, we illustrate and 
quantify the rationale behind for error suppression and 
its promising exploitation potential. Our approach is also 
consistent with the fact that none of the 9 13 experiments 
have so far provided detailed publications on the non¬ 
trivial reactor systematics analyses, including the quan¬ 
titative justification of all assumptions used so far. In 
the context of the $13 experiments, a few references exist 
DC [161119] and DB[20] (nothing yet available for RENO) 
illustrating a fraction of their reactor flux studies, but not 
dealing with the inter-reactor error correlation needed for 
the exploitation of SF(correction). DC is, however, final¬ 
ising a dedicated publication m on their reactor flux 
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systematics quoted so far. Thus, reactor flux error crit¬ 
ical subject, despite its major impact to the final pre¬ 
cision on 0 i 3 , remains unfortunately somewhat obscure 
in today’s literature. As time goes, the statistical error 
of reactor experiments is no longer dominant, the treat¬ 
ment of the systematics should be clearly laid well in 
advance to maximise the stability and reliability of the 
013 measurement, whose impact has critical implication 
transcending reactor neutrino results, affecting, for ex¬ 
ample, current searches for neutrino CP-violation. 


CONCLUSIONS 

This publication provides reactor neutrino experiments 
with a coherent treatment for reactor flux systematic un¬ 
certainty for multi-detector experiments in multi-reactor 
contributions. We started with careful treatment of the 
single detector in a single reactor site scenario, being 
the most relevant case for most reactor experiments be¬ 
yond 0i 3 . We have demonstrated that the challeng¬ 
ing reactor flux systematic do not trivially cancel by 
the adoption of that multi-detector experiments. How¬ 
ever, we have identified several means for error suppres¬ 
sion in the context of Double Chooz, Daya Bay and 
RENO experiments, using simplified scenarios to max¬ 
imise relative comparison. We computed an integral er¬ 
ror suppression fraction (SF), which can be broken down 
into three components, defined as SF(total)=SF(iso- 
flux)xSF(A' fl )xSF (correlation), where SF (Nr) sup¬ 
presses the uncorrelated error of identical reactors, 
SF (iso-flux) suppresses the total error if the site geome¬ 
try meets fully, or partially, the iso-flux condition and 
SF(correlation) suppresses the error if the reactor er¬ 
rors are fully correlated. SF(iso-flux) and SF (correlation) 
could lead to total cancellation of the flux error. How¬ 
ever, only SF (Nr) has been exploited to improve the 0i3 
precision, although total cancellation is impossible. 

This publication deals in detail on the calculation for 
SF (iso-flux), thus paving the ground for its exploitation, 
yielding two important observations. First, DC, once in 
its near-(-far configuration, is the only experiment ex¬ 
pected to benefit from a negligible reactor flux error, 
thanks to the ~90% iso-flux error suppression. Second, 
Daya Bay and RENO could also benefit from their partial 
iso-flux, thus yielding up to ~50% flux error suppression. 
Thus, this publication embodies a major improvement in 
the global precision of 0 i 3 by improving the precision of 
all experiments measuring it, including current results. 
Finally, we have highlighted the potential for a mecha¬ 
nism, currently neglected, for error suppression relying 
on further reactor error correlation insight, characterised 
by the SF(correlation) term, further improving the pre¬ 
cision of all multi-detector experiments. 
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Appendix: Error Propagation for Partially 
Correlated Uncertainties 

In order to deliver the correlation coefficient for par¬ 
tially correlated uncertainties, we generalise the approach 
presented in [18]. Keeping the same notations, let us 
consider a measurement F which depends on two vari¬ 
ables x and y having the absolute uncertainties: Aa; and 
Ay. We split one of them, for example Ax, into its 
two components Ax u and Ax c representing respectively 
the totally uncorrelated and totally correlated both rel¬ 
ative to Ay. Since the general error propagation for¬ 
mula is symmetric and the correlation is mutual, one 
can split whether Ax or Ay. The fraction between the 
correlated and uncorrelated components is given charac¬ 
terised by a constant a such that Ax u = aAx c . Since 
(Aa:) 2 = (Ax u ) 2 + (Ax c ) 2 , the uncorrelated and corre¬ 
lated components are expressed as 

. ,, aAx , . _ Aa; .. 

Aa; = , and Aa; = , (A26) 

v / ?TT v^TT 

In order to estimate the total uncertainty 

(A F) 2 = (A F u f + (A F c ) 2 (A27) 

we shall compute each of its components. The uncorre¬ 
lated component is 


(A F u ) 2 



2 

(Aa; 11 ) 2 


(Aa 


(A28) 


while the correlated component is 


(AF C ) 2 = 


£)W + (£)W + 
^(£)W+(f)W + 


2 f8F\ fdF\ 

WTT UJ Al U 1A " 


(A29) 


Replacing Eq. |A28| and Eq. |A29| in Eq. |A27| we obtain 


< af ) 2 = ('£)W + (f) 2 (^> 2+ 

ysMIMf)*' (A30) 


Comparing Eq. |A30| to the general law of error propaga¬ 
tion 


< aF > 2 ^ (S) 2 ^) 2 + (f ) 2(A!,)2 + 

+2k (%) Av (A31) 


we define the correlation error coefficient as 


k = 


Va 2 + 1 


(A32) 


Replacing now the a parameter, we can rewrite the pre¬ 
vious equation as 


k = 


1 + 


Ax c —Ax u 

Ax c -\-Ax u 


( Ax c — Ax u \ 
\ Ax c +Ax u J 


(A33) 


Eq. |A33| remains the same for the relative uncertainty 


k = 


1 + 


s c -s u 

<5 c +<5“ 




(A34) 


a 2 + 1 


OF 

dx 
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where S u,c stands for Ax u,c /x, respectively. 




































